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CTs 1 Abstract 

o 



Interactions of virtual photons are analyzed in terms of photon structure. It is argued 
that the concept of parton distribution functions is phenomenologically very useful 
even for highly virtual photons involved in hard collisions. This claim is illustrated on 
Q | leading order expressions for F^(x, P 2 , Q 2 ) and effective parton distribution function 

f-^. ' D c g(x, P 2 ,Q 2 ) relevant for jet production, as well as within the next -to-leading order 

CN . QCD calculations of jet cross-sections in electron-proton collisions. 

^ ■ 1 Introduction 

Parton distribution functions (PDF) are, together with the colour coupling a s , the basic ingredients 
of perturbative QCD calculations. It is worth emphasizing that in quantum field theory it is difficult 
Q\ ■ to distinguish effects of the "structure" from those of "interactions" . Within the Standard Model 

(SM) it makes good sense to distinguish fundamental particles, which correspond to fields in its 
^ lagrangian jCsm (leptons, quarks and gauge bosons) from composite particles, which appear in the 

mass spectrum but have no corresponding field in £sm- For the latter the use of PDF to describe 



Oh 
D 



their "structure" appears natural, but the concept of PDF turns out to be phenomenologically 
useful also for some fundamental particles, in particular the photon. PDF are indispensable for the 



real photon due to strong interactions between the qq pair to which it couples electromagnetically. 
For massless quarks this coupling leads to singularities, which must be absorbed in PDF of the 
photon, similarly as in the case of hadrons. Although PDF of the real photon satisfy inhomogeneous 
evolution equations, their physical meaning remains basically the same as for hadrons. 

For nonzero photon virtualities there is no true singularity associated with the coupling 7* — > qq 
and one therefore expects that for sufficiently virtual photon its interactions should be calculable 
perturbatively, with no need to introduce PDF. The main aim of this paper is to advocate the use 
of PDF also for virtual photons involved in hard collisions. 

Throughout this paper we shall stay within the conventional approach to evolution equations for 
PDF of the photon. The reformulation of the whole framework for the description of hard collisions 
involving photons in the initial state, proposed recently by one of us Q, affects the analysis of dijet 
production in photon-proton collisions, which is the main subject of this paper, only marginally. 
This stands in sharp contrast to QCD analysis of F7 which is affected by this reformulation quite 
significantly. 

The paper is organized as follows. In Section 2 the notation and basic facts concerning the 
evolution equations for PDF of the real photon and the properties of their solutions are recalled. In 
Section 3 the physical content of PDF of the virtual photon is analyzed and its phenomenological 
relevance illustrated within the LO QCD. This is further underlined in Section 4 by detailed analysis 
of NLO QCD calculations of dijet cross-sections in ep collisions, followed by the summary and 
conclusions in Section 5. 
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2 PDF of the real photon 



In QCD the coupling of quarks and gluons is characterized by the renormalized colour coupling 
("couplant" for short) a s (p), depending on the renormalization scale \i and satisfying the equation 

^ ee flco*)) = -^o?M - ^o?M + -.., (1) 

where, in QCD with rif massless quark flavours, the first two coefficients, (3q = 11 — 2ny/3 and 
Pi = 102 — 38n//3, are unique, while all the higher order ones are ambiguous. As we shall stay 
in this paper within the NLO QCD, only the first two, unique, terms in (|l]) will be taken into 
account. However, even for a given r.h.s. of ([!]) its solution a s (fx) is not a unique function of /i, 
because there is an infinite number of solutions of ([!]), differing by the initial condition. This so 
called renormalization scheme (RS) ambiguity [] can be parameterized in a number of ways. One of 
them makes use of the fact that in the process of renormalization another dimensional parameter, 
denoted usually A, inevitably appears in the theory. This parameter depends on the RS and at the 
NLO even fully specifies it. For instance, a s (/i) in the familiar MS and MS RS are two solutions 
of the same equation (|]), associated with different Ars- At the NLO the variation of both the 
renormalization scale [i and the renormalization scheme RS={Ars} is redundant. It suffices to fix 
one of them and vary the other, but we stick to the common practice of considering both of them 
as free parameters. In this paper we shall work in the standard MS RS of the couplant. 

The "dressed" PDF result from the resummation of multiple parton collinear emission off 
the corresponding "bare" parton distributions. As a result of this resummation PDF acquire 
dependence on the factorization scale M. This scale defines the upper limit on some measure t of 
the off-shellness of partons included in the definition of D(x, M) 

rM 2 

Di(x,M)= / dtdi(x,t), i = q,q,G, (2) 

where the unintegrated PDF di(x,t) describe distributions of partons with the momentum fraction 
x and fixed off-shellness t. Parton virtuality r =| p 2 — m 2 | or transverse mass m|i = p\ + m 2 , are 
two standard choices of such a measure. Factorization scale dependence of PDF of the photon is 
determined by the system of coupled inhomogeneous evolution equations 



dS(x, M) 
dlnM 2 
dG{x,M) 



8nk q + P qq ®Y, + P qG ® G, (3) 

dlnM 2 " k G + P Gq ®X + P GG ®G, (4) 
dq NS (x,M) 

dlnM2 = oviskq + Pns ® <?NS, (5) 



where 



n; n; 

E(x,M) ee ^q+(x,M)=J2k&M)+q i (x,M)}, (6) 

i=l i=l 
n f 

q NS (x,M) ee J2(^-( e2 ))(li(x,M)+q i (x,M)), (7) 
i=l 



x In higher orders this ambiguity includes also the arbitrariness of the coefficients > 2. 

2 In the following the adjective "dressed" will be dropped, and if not stated otherwise, all PDF will be understood 
to pertain to the photon. 
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^NS = 6n / ((e 4 )-(e 2 ) 2 ), <5 s = 6n / (e 2 ). 
To order a the splitting functions Pij and ki are given as power expansions in a s (M): 

' k?){x) + sum k iD {x) + ( ^m\v{ X) + . . 



k q {x,M) = — 



k G (x,M) 
Pij(x,M) 



2tt 

a 
2^ 



2tt 



2tt 
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a,{M) 
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) 4 2) (*) + 



°^P, ( ° ) (x) + 



2vr « 



2tt 



(8) 

0) 
(10) 
(11) 



where the leading order splitting functions k q 



x 



+ (1 



and P^{x) are unique, while 



k%\Pjf, ,j > 1 depend on the choice of the factorization scheme (FS). 
i as evolutior 
3c,' A?,? 



all higher order ones kq\ltQ , 

The equations (||-||) can be rewritten as evolution equations for qi(x, M),q i (x, M) and G(x,M) 
with inhomegenous splitting functions kq^ = 3e 2 k q °\ The photon structure function F^ix^Q 2 
measured in deep inelastic scattering of electrons on photons is given as a sum of convolutions 



x 



a 



g N s(M) ® C q (Q/M) + — 5 NS C 7 + 

Z7T 



(12) 



a 



(e 2 )S(M) ® C,(Q/M) + — (e 2 )feC 7 + < e ^G(M) ® C G (Q/M) (13) 



2vr 



of photonic PDF and coefficient functions C q (x), C G (x), C 7 (x) admitting perturbative expansions 



C q (x,Q/M) 
C G {x,Q/M) 
C y (x,Q/M) 



S(l-x) + ^M C W(x,Q/M) + --, 
^Cj>\x,Q/M) + ..., 
CV>{x, Q/M) + %&C&)(z, Q/M) + ■■-. 



(14) 
(15) 
(16) 



The renormalization scale fi, used as argument of ct s (ii) in (14-16) is in principle independent of 
the factorization scale M. Note that despite the presence of as argument of a s (/j,) in (14-ilf), 
the coefficient functions C q , C G and C 7 , summed to all orders of a s , are actually independent of it, 
because the /i-dependence of the expansion parameter a s (fi) is cancelled by explicit dependence 

(i) (i) (i) 

of Cq . C 1 q , C 7 , 2 ^ 2 on On the other hand, PDF as well as the coefficient functions C q ,C G 
and C 7 do depend on both the factorization scale M and factorization scheme, but in such a 
correlated manner that physical quantities, like F% , are independent of both M and the FS, provided 
expansions (PHTll) and (|l4]-p^) are taken to all orders in a s (M) and a s (fi). In practical calculations 
based on truncated forms of (p|-[il|) and (14-lq) this invariance is, however, lost and the choice of 
both M and FS makes numerical difference even for physical quantities. The expressions for 
C^ , Cq given in || are usually claimed to correspond to "MS factorization scheme" . As argued 
in [H], this denomination is, however, incomplete. The adjective "MS" concerns exclusively the 
choice of the RS of the couplant a s and has nothing to do with the choice of the splitting functions 
pjp . The choices of the renormalization scheme of the couplant a s and of the factorization scheme 
of PDF are two completely independent decisions, concerning two different and in general unrelated 
redefinition procedures. Both are necessary in order to specify uniquely the results of fixed order 
perturbative calculations, but we may combine any choice of the RS of the couplant with any 
choice of the FS of PDF. The coefficient functions C q ,C G ,Cy depend on both of them, whereas 
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G PL (x,M 2 ) 

X,T 




Figure 1: Diagrams denning the pointlike parts of nonsinglet quark and gluon distribution functions. 
The resummation involves integration over parton virtualities r < M 2 . 

the splitting functions depend only on the latter. The results given in Q correspond to MS RS of 
the couplant but to the "minimal subtraction" FS of PDF 0. It is therefore more appropriate to 
call this full specification of the renormalization and factorization schemes as "MS + MS scheme" . 
Although the phenomenological relevance of treating [i and M as independent parameters has been 
demonstrated [||, we shall follow the usual practice of setting \x = M. 

2.1 Pointlike solutions and their properties 

The general solution of the evolution equations (||-[|) can be written as the sum of a particular 
solution of the full inhomogeneous equation and the general solution of the corresponding homo- 
geneous one, called hadronic Q part. A subset of the solutions of full evolution equations resulting 
from the resummation of series of diagrams like those in Fig. [I], which start with the pointlike 
purely QED vertex 7 — * qq, are called pointlike (PL) solutions. In writing down the expression 
for the resummation of diagrams in Fig. [I] there is a freedom in specifying some sort of boundary 
condition. It is common to work within a subset of pointlike solutions specified by the value of the 
scale Mq at which they vanish. In general, we can thus write (D = q, q, G) 

D{x, M 2 ) = D PL (x, M 2 ) + D HAD (x, M 2 ). (17) 

Due to the fact that there is an infinite number of pointlike solutions q PL (x,M 2 ), the separation 
of quark and gluon distribution functions into their pointlike and hadronic parts is, however, am- 
biguous and therefore these concepts have separately no physical meaning. In we discussed 
numerical aspects of this ambiguity for the Schuler-Sjostrand sets of parameterizations JjJ. 

To see the most important feature of the pointlike part of quark distribution functions that 
will be crucial for the following analysis, let us consider in detail the case of nonsinglet quark 
distribution function q^s(x,M), which is explicitly defined via the series 

2vr 3 J M 2 t J x y qq \yj Jm§ n 2vr 2vr 3 J M 2 t 2 

3 See Section 2.6 of 0, in particular eq. (2.31), for discussion of this point. 

4 Sometimes also called "VDM part" because it is usually modelled by PDF of vector mesons. 
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Figure 2: Comparison, left in the nonsinglet and right in the singlet channels, of the functions 
xki(x) = xSikq ,i = NS,S with the function (f2~l|) and its analogue in the singlet channel. 



p(0) 
99 



1 dw p(Q) (?L\ f M * dn a s (ri) f Tl dr 2 a s (r 2 ) a j ) 
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where fc^g(x) = ^NS^g ^)- I n terms of moments defined as 



k ( >(w) 



M 2 r 3 



+ •••, (18) 



f(n) = / x n f(x)dx 
Jo 



(19) 



this series can be resummed in a closed form 

4tt 



q^{n,M ,M) 



a s {M) 



1 



a 



l-2P^ ) (n)//3 



a s (M )J 



ONs(n), 



(20) 



where 



a 



27T/3oi_ 2 P (? ( g 0) (n)//3o 



(21) 



It is straightforward to show that (|T^) or, equivalently, (^0|), satisfy the evolution equation (||) 



with the splitting functions k q and Pj, including the first terms kq and Pgg'' only. 

Transforming (^Tj) to the x -space by means of inverse Mellin transformation we get ons^) 
shown in Fig. 2. The resummation softens the x— dependence of clns( x ) with respect to the first 
term in fll8|), proportional to kj<fs(x), but does not change the logarithmic dependence of c/ns on 
M. In the nonsinglet channel the effects of gluon radiation on q^g are significant for x > 0.6 but 
small elsewhere, whereas in the singlet channel such effects are marked also for x < 0.5 and lead to 
a steep rise of xq^g at very small x. As emphasized long time ago by authors of [|| the logarithmic 
dependence of g^g on InM has nothing to do with QCD and results exclusively from integration 
over the transverse momenta (virtualities) of quarks coming from the basic QED 7* — > qq splitting. 
For M/Mq — > 00 the second term in brackets of (|20| ) vanishes and therefore all pointlike solutions 
share the same large M behaviour 



,(o) 



q^{x,M ,M) 



47T 



ONS(Z) = <?NS ( X ' M ) a ln 



a,{M) 



A 2 



(22) 
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Figure 3: The u and c quark and gluon distribution functions of the real photon for SaSID (upper 
solid curves) and SaS2D (upper dotted curves) parameterizations at M 2 = 100 GeV 2 . The VDM 
and pointlike parts of both parameterizations are plotted separately, the latter as solid and dotted 
curves peaking at large x. For quarks the splitting terms (|23|) corresponding to SaSID and SaS2D 



(upper and lower dashed curves) are overlayed to show the effects of resummation (18). 



defining the so called asymptotic pointlike solution q^(x,M) The fact that for the asymp- 

totic pointlike solution ( p2|) a s (M) appears in the denominator of ( p2|) has been the source of 
misleading claims that q(x,M) = 0(a/a s ). In fact, as argued in detail in ||], q(x,M) = 0(a) 
as suggested by explicit construction in (^). We shall return to this point in Section 4.2 when 
discussing the factorization scale invariance of finite order approximations to dijet cross-sections. 

The arbitrariness in the choice of Mq reflects the fact that as Mq increases, less of the gluon 
radiation effects is included in the resummation (|20| ) defining the pointlike part of quark distribution 
function g PL but included in hadronic one. The latter is usually modelled by the VDM ansatz and 
will therefore be called "VDM" in the following. As we shall see, the hadronic and pointlike parts 
have very different behaviour as functions of x and M. 



2.2 Properties of Schuler— Sjostrand parameterizations 

Practical aspects of the ambiguity in separating PDF into their VDM and pointlike parts can be 
illustrated on the properties of SaSID and SaS2D parameterizations |Q] 0. What makes the SaS 
approach particularly useful for our discussion is the fact that it provides separate parameterizations 
of the VDM and pointlike parts [] of both quark and gluon distributions. 

In Fig. H distribution functions xu(x, M 2 ), xc(x,M 2 ) and xG(x,M 2 ) as given by SaSID and 
SaS2D parameterizations for M 2 = 100 GeV 2 are compared. To see how much the resummation 
of gluon radiation modifies the first term in ( |l8|) we also plot the corresponding splitting terms 

M 2 ,M 2 ) = |-3e 2 (x 2 + (1 - x) 2 ) In (23) 

In Fig. |i] the scale dependence of VDM and pointlike parts of u and c quark and gluon distribution 
functions is displayed. In the upper six plots we compare them as a function of x at M 2 = 
25, 100, 1000 GeV 2 , while in the lower three plots the same distributions are rescaled by the factor 

5 The properties of SaSIM and SaS2M parameterizations are similar. 
6 Called for short "pointlike quarks" and "pointlike gluons" in the following. 
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Figure 4: Factorization scale dependence of parton distributions functions u(x,M),c(x,M) and 
G(x, M) of the real photon. Dashed and solid curves correspond, in the order indicated by the 
arrows, to pointlike and VDM parts of these distributions at M 2 = 25, 100 and 1000 GeV 2 . In the 
lower part SaSID quark distribution functions xu(x,M 2 ,P 2 ) rescaled by ln(M 2 /MQ) for the real 
photon (left) and by ln(M 2 /P 2 ) for the virtual one, are plotted and compared to the predictions 
of the splitting term (E9). 
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Figure 5: F^{x, Q 2 ) as given by SaSID (solid curves) and SaSD2 (dotted curves) parameterizations. 
The full results are given by the upper, the pointlike and VDM contributions parts by two lower 
curves. The dashed curves describe the contributions of the splitting term (23). 
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Figure 6: D e g(x, M 2 ) as given by SaSID and SaSD2 parameterizations. Solid curves show the full 
results, dashed ones the VDM contributions. The pointlike parts are separated into the contribu- 
tions of pointlike quarks (dotted curves) and pointlike gluons (dash-dotted). Thick dashed curves 
correspond to the splitting term ( |23| ) . 
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ln(M 2 /Mq). Dividing out this dominant scale dependence allows us to compare the results directly 
to curves in Fig. ^. Finally, in Figs. |B| and |6| SaS predictions for two physical quantities, F^ and 
effective parton distribution function -D e ff > relevant for jet production in 7p and 77 collisions, 

Fl(x,Q 2 ) = ^2^ g (x,Q 2 ), (24) 
D cS (x,M 2 ) = ^( % (x,M 2 )+^(x,M 2 )) + ^G(x,M 2 ) (25) 

2=1 

are displayed and compared to results corresponding to the splitting term (23). Figures |3|-|6] illus- 
trate several important properties of PDF of the real photon: 

• There is a huge difference between the importance of the VDM components of light quark 
and gluon distribution functions in SaSID and SaS2D parameterizations: while for SaS2D 
the VDM component is dominant up to x = 0.75, for SaSID the pointlike one takes over 
above x = 0.1! 

• Factorization scale dependence of VDM and pointlike parts differ substantially. While the 
former exhibits the pattern of scaling violations typical for hadrons, the latter rises, for quarks 
as well as gluons, with M for all x. For pointlike gluons this holds despite the fact that 
G PL (x,M 2 ) satisfies at the LO standard homogeneous evolution equation and is due to the 
fact that the evolution of G PL (x, M 2 ) is driven by the corresponding increase of E PL (x,M 2 ). 

• The approach of g PL (x,M 2 ) to the asymptotic pointlike solution (22) is slow []. For experi- 
mentally relevant values of M and except for x close to 1 the effects of multiple gluon emission 
on pointlike quarks are thus small. 

• As the factorization scale M increases the VDM parts of both quark and gluon distribution 
functions decrease relative to the pointlike ones, except in the region of very small x. 

• Despite huge differences between SaSID and SaS2D parameterizations in the decomposition 
of quark and gluon distributions into their VDM and pointlike parts, their predictions for 
physical observables F7 and D e s are much closer. 



The most prominent effects of multiple parton emission on physical quantities appear to be 
2 

sections, approximately described by D e g(x, M 2 ). 



the behaviour of F?(x,Q 2 ) at large x and the contribution of pointlike gluons to jet cross- 



3 PDF of the virtual photon 

For the virtual photon the initial state singularity due to the splitting 7* — > qq is shielded off by 
the nonzero initial photon virtuality P 2 and therefore in principle the concept of PDF does not 
have to be introduced. In practice this requires, roughly, P 2 > 1 GeV 2 , where perturbative QCD 
becomes applicable. Nevertheless, PDF still turn out to be phenomenologically very useful because 

• their pointlike parts include the resummation of parts of higher order QCD corrections, 

• the hadronic parts, though decreasing rapidly with P 2 , are still vital at very small x. 

7 The fact that for the c-quark the SaS2D parameterization is close to the curve corresponding to g spllt (:r, M§, M 2 ), 
whereas SaSID lies substantially below, reflects the fact that Mo > m c for SaS2D, but Mo < m c for SaSID 
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Both of these aspects define the "nontrivial" structure of the virtual photon in the sense that they 
are not included in the splitting term ( |23| ) and thus are not part of existing NLO unsubtracted 
direct photon calculations. One might argue that the calculable effects of resummation should 
be considered as "interaction" rather than "structure", but their uniqueness makes it natural to 
describe them in terms of PDF. 

3.1 Equivalent photon approximation 

All the present knowledge of the structure of the photon comes from experiments at the ep and 
e + e~ colliders, where the incoming leptons act as sources of transverse and longitudinal virtual 
photons. To order a their respective unintegrated fluxes are given as 



The transverse and longitudinal fluxes thus coincide at y = 0, while at y = 1, f? vanishes. The 
1/.P 2 dependence of the first terms in (26-27) results from the fact that in both cases the vertex 



where photon is emitted is proportional to P 2 . This is due to helicity conservation for the transverse 
photon and gauge invariance for the longitudinal one. The term proportional to m 2 /P 4: in ( p6| ) 
results from the fact that the helicity conservation at the e7e vertex is violated by terms proportional 
to electron mass. No such violation is permitted in the case of gauge invariance, hence the absence 
of such term in ([27]). Note that while for P 2 S> m 2 the second term in ( f26| ) is negligible with 
respect to the leading l/P 2 one, close to -P min = m 2 y 2 /(l — y) their ratio is finite and approaches 
2(l-y)/(l + (l-y) 2 ). 

3.2 Lessons from QED 

The definition and evaluation of quark distribution functions of the virtual photon in pure QED 
serves as a useful guide to parton model predictions of virtuality dependence of the pointlike part 
of quark distribution functions of the virtual photon. In pure QED and to order a the probability 
of finding inside the photon of virtuality P 2 a quark with mass m q , electric charge e g , momentum 
fraction x and virtuality r = m 2 — k 2 < M 2 is defined as (k denotes its four-momentum) 



q Q En(x,m 2 q ,P 2 ,M 2 )= (£3e 2 ) / _^_i^dr, (28) 



_ (^A [ M2 w ( x > m l p 

where the function W(x,m 2 ,P 2 ,M 2 ) can in general be written as 

W{x,m 2 ,P 2 ,M 2 ) = f{x)Y^-+g{x)m 2 q + h{x)P 2 + ••• 



f(x)r + ( g(x) - ^4 ) m 2 + (h{x) - xf(x)) P 2 + • • • . (29) 



m 

1 — x 

In the collinear kinematics, which is relevant for finding the lower limit on r, the values of m q , x, r 
and pt are related to initial photon virtuality P 2 as follows 

„ m 2 + p% „ m 2 

T = xP 2 + q PT , => T mm = xP 2 + q -. (30) 

1 — x 1 — x 
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The functions f(x),g(x) and h(x), which determine the terms singular at small r, are unique 
functions that have a clear parton model interpretation: so long as r < M 2 eq. fl2q) describes the 
flux of quarks that are almost collinear with the incoming photon and "live" long with respect to 
1/M. On the other hand, the terms indicated in (|29| ) by dots are of the type T k+1 /s k , k > 1, which 
upon insertion into ( |28| ) yield contributions that are not singular at r = and therefore do not 
admit simple parton model interpretation. In principle we can include in the definition ( |28| ) even 
part of these nonpartonic contributions, but we prefer not do that. Substituting (p9| ) into (pg|) and 
performing the integration gives, in units of 3e 2 a/27r, 



q Q ED(x,m 2 q ,P 2 ,M 2 ) = f(x) In 



M 2 



+ 



Six) + 



g(x)m 2 + h(x)P 2 



M 2 



(31) 



In practical applications the factorization scale M is identified with some kinematical variable 
characterizing hardness of the collision, like sfQ 2 in DIS or Ej, in jet production. For r mm <C M 2 
the expression simplifies to 



q Q ED(x,m 2 q ,P 2 ,M 2 ) =f(x)]n( 



M 1 



xP 2 + m 2 /(l — x) 



g(x)vn 2 + h(x)P 2 
xP 2 + m 2 /(l — x) 



(32) 



For x(l — x)P 2 3> m 2 this expression reduces further to 



q Q Ev(x,0,P 2 ,M 2 ) =/(x)ln 



M 2 

xT 2 



f(x) + 



h(x) 



Provided m 2 ^ (|32]) has a finite limit for P 2 — » 0, corresponding to the real photon 



(33) 



q QE D(x,m 2 q ,0,M 2 ) = f(x)ln 



M 2 (l 



in,, 



Six) +g(x)(l - x). 



(34) 



As in the case of the photon fluxes (|2^-p7|) the logarithmic term, dominant for large M 2 , as well 
as the "constant" terms, proportional to f(x),g(x) and h(x), come entirely from the integration 
region close to T mm and are therefore unique. At r = r mm both types of the singular terms, i.e. 
1/r or 1 /t 2 , are of the same order but the faster fall-off of the 1/r 2 terms implies that for large M 2 
the integral over r, which gives (|3l|), is dominated by the weaker singularity 1/r. In other words, 
while the logarithmic term is dominant at large M 2 , the constant terms resulting from nonzero m 2 
and P 2 come from the kinematical configurations which are even more collinear, and thus more 
partonic, than those giving the logarithmic term. The analysis of the vertex 7* — > qq in collinear 



kinematics yields [12| 



f T (x) = x 2 + (1 -x) 2 , g T {x) 
Sl(x) = 0, 



1 



-, hr(x) = Q, 



1 — x ' 

g L ix) = 0, h L {x)=Ax 2 {l-x). 



(35) 



The expressions (|3l|-p2;) exhibit explicitly the smooth transition between quark distribution func- 
tions of the virtual and real photon. This transition is governed by the ratio P 2 /m 2 , which un- 
derlines why in QED fermion masses are vital. On the other hand, as P 2 (or more precisely 
r min ) increases toward the factorization scale M 2 , the above expressions for the quark distribution 
functions of the virtual photon vanish. This property holds not only for the logarithmic term 
but also for the "constant" terms and has a clear intuitive content: virtual photon with lifetime 
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1 /P <C 1 /M does not contain partons living long enough to take part in the collision characterized 
by the interaction time 1/M. 

For virtual photon and x(l-x)P 2 > m 2 the coefficient functions C^(x, P 2 ,Q/M),C^(x, P 2 , Q/M) 
for transverse and longitudinal target photon polarizations are given as |13| 



C^(x,P 2 ,l) 
C^P 2 ,!) 



(x 2 + (1 - xf) In + 8x(l - x) - 2 
x z 



: 4x(l — X), 

whereas for the real photon, i.e. for P 2 = 



c$M,i) 

C^°i(x,0,l) 



(x 2 + (1 - x) 2 )ln 



1 — x 



+ 8x(l -x) - 1 



(36) 
(37) 

(38) 
(39) 



The origins of the nonlogarithmic parts of in (|36-38) can then be identified as follows: 
-l + 8x(l-x) = -2 + 8x(l-x)+ 1 

> v ' 

for massless quark 9t 

(x)(l-x) 

= -l + 6x(l-x)-(x 2 + (l-x) 2 )+ 1 

g T (x)(l-x) 



(40) 



nonpartonic part 



fr(x) 



The nonpartonic part itself can be separated into two pieces, coming from the interaction of the 
transverse target photon with transverse and longitudinal probing one 



1 + 6x(l — x) 



-1 + 2x(l - x) -4x(l - x) 

S V - 

from (Ttt from o"lt 



(41) 



Except for a brief comment in the next subsection, we shall consider throughout the rest of this 
paper the contributions of the transverse polarization of the target virtual photon only. The 
importance of including in analyses of hard collisions of virtual photons the contributions of 7^ will 



be discussed in detail in separate publication [16]. 



3.3 What is measured in DIS on virtual photons? 

In experiments at e + e~ colliders the structure of the photon has been investigated via standard 
DIS on the photon with small but nonzero virtuality P 2 . The resulting data were used in [14] to 
determine PDF of the virtual photon. In these analyses was taken in the form 



Cf(x,P 2 ,l) = 3 



1 



(x z + (1 - xf) In + 6x(l - x) - 2 



(42) 



which, however, does not correspond to the structure function that is actually measured in e + e 
collisions, but to the following combination 



F^(x,P 2 ,Q 2 ) = Fl T (x,P 2 ,Q 2 ) - -Fl L (x,P 2 ,Q 2 



(43) 



of structure functions corresponding to transverse and longitudinal polarizations of the target pho- 
ton. This combination results after averaging over the target photon polarizations by means of 
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contraction with the tensor —g^ u /2. The expression —2 + 6x(l — x) follows also directly from the 
definition ([4S|) and considerations of the previous Section: 

- 2 + 6x(l - x) = -2 + 8s(l - s) - 2x(l - x) . (44) 

V v 

Iron 7 T from 7l/2 



Because the fluxes (26|-27) of transverse and longitudinal photons are different functions of y, 
any complete analysis of experimental data in terms of the structure functions F^ T {x,P 2 ,Q 2 ) 
and F^ L (x, P 2 ,Q 2 ) at fixed x,P 2 ,Q 2 requires combining data for different y. This is in prin- 
ciple possible, but experimentally difficult to accomplish. The situation is simpler at small y, 
where f^(y,P 2 ) = f2(y,P 2 ) = f 1 {y,P 2 ), and the data therefore correspond to the convolution 
of / 7 (y, P 2 ) with the sum F^ T + F2 L . ^he nonlogarithmic term in corresponding to this 



combination is, however, not —2 + 6x(l — x), as used in [14]], but —2 + 12x(l — x), the sum of non- 
logarithmic terms corresponding to transverse and longitudinal photons. Numerically the difference 
between these two expressions is quite sizable. 

Very recently, the GRS group [15| has changed their approach to the treatment of the target 
photon polarizations and argued in favor of neglecting the contribution of 7^ and using even for 

the virtual photon the same form of as for the real one. We disagree with their arguments, 
but leave the discussion of this point to future paper [If 



3.4 Virtuality dependent PDF 

In realistic QCD the nonperturbative effects connected with the confinement, rather than current 
or constituent quark masses, are expected to determine the long-range structure of the photon 
and hence also the transition from the virtual photon to the real one. For instance, within the 
SaS parameterizations the role of quark masses is taken over by vector meson masses for the VDM 
components and by the initial Mq for the pointlike ones. As in the case of the real photon, we 
recall basic features of SaS parameterizations of the virtual photon, illustrated in Figs. |7|-|l0|. 

• With increasing P 2 the relative importance of VDM parts of both quark and gluon distribution 
functions with respect to the corresponding pointlike ones decreases rapidly. For M 2 > 25 
GeV 2 the VDM parts of both SaSID and SaS2D parameterizations become negligible for 
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Figure 8: The same as in Fig. |] but for virtual photon with P 2 = 3 GeV 2 . 
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Figure 9: Virtuality dependence of photonic PDF at M 2 = 100 GeV 2 . Dashed curves correspond 
to the pointlike and solid ones to the VDM parts, from above for P 2 = 0, 1, 3 GeV 2 . 
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Figure 10: The same as in Fig. |B| but for virtual photon with P 2 = 3 GeV 2 . 
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Figure 11: D e g(x, P 2 , M 2 ) as a function of x for M 2 
Notation as in Fig. |6|. In the splitting term Mq = P 2 . 



25,100,1000 GeV 2 and P 2 = 3 GeV 2 . 
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P 2 > 2 GeV 2 , except in the region of very small x < 0.01. Consequently, also the ambiguity 
in the separation (Il7|) is practically negligible in this region. 



The general pattern of scaling violations, illustrated in Fig. [8], remains the same as for the real 
photon, but there is a subtle difference, best visible when comparing in Fig. |I] the rescaled 
PDF for P 2 = with those at P 2 = 1,3 GeV 2 . While for P 2 = 0, increasing M 2 softens the 
spectrum towards the asymptotic pointlike form (^), for P 2 > 1 GeV 2 quark distribution 
functions increase with M for practically all values of x. Moreover, while for P 2 = the 
splitting term intersects the SaSID curves qualitatively as in Fig. for P 2 > 1 GeV 2 it is 
above them for all x . These properties reflect the fact that SaS parameterizations of PDF of 
the virtual photon do not satisfy the same evolution equations as PDF of the real one. This 
difference is formally of power correction type and thus legitimate within the leading twist 
approximation we are working in, but numerically nonnegligible. 

As shown in Fig. ^ both VDM and pointlike parts decrease with increasing P 2 , but the VDM 
parts drop much faster than the pointlike ones. 



The implications of these properties for physical quantities F^x, P 2 , Q 2 ) and D e g (x, P 2 , M 2 ) are 
illustrated by Figs. 11 and 1C. For F^ix, P 2 ,M 2 ) the nontrivial aspects of virtual photon structure 



are confined mainly to the region of large x, where they reduce the predictions based on the 
splitting term (23). Obviously, only the pointlike quarks are relevant for this effect. The effects on 



D c ff(x,P , M ) worth emphasizing are the following: 

• For all scales M the contribution of the splitting term is above the one from pointlike quarks, 
the gap increasing with increasing P 2 and decreasing M 2 . 

• Pointlike quarks dominate D c g (x, P 2 , M 2 ) at large x, whereas for x < 0.3, most of the point- 
like contribution comes from the pointlike gluons. In particular, the excess of the pointlike 
contributions to D c g over the contribution of the splitting term, observed at x < 0.5, comes 
almost entirely from the pointlike gluons! 



• For x > 0.6 the full results are clearly below those given by the splitting term (|23| ) with 
Mg = P 2 . In this region one therefore expects the sum of subtracted direct and resolved 
contributions to jet cross-sections to be smaller than the results of unsubtracted direct cal- 
culations. 

So far no restrictions were imposed on transverse energies and pseudorapidities of jets produced 
in 7*p collisions. However, as discussed in more detail in Section 4.3.4, for jet transverse energies 
Ft ^ 5 GeV, hadronization corrections become intolerably large and model dependent in the 
region r/ < —2.5. On the experimental side the problems with the reconstruction of jets with 
Ft ^ 5 GeV in the forward region restrict the accessible region to ry < 0. The cuts imposed on 
pseudorapidities r]^\i = 1,2 and transverse energies Ej) of two final state partons ^ influence 
strongly the corresponding distribution in the variable x^ 

" s ~ (45) 

used in analyses of dijet data. At the LO and for massless partons x 7 defined in (fl5|) coincides with 
the conventional fraction x appearing as an argument of photonic PDF. MC simulations show that 

8 All quantities correspond to 7*p cms. 

9 In realistic QCD analyses of two jets with highest and second highest Et- Jets with highest and second highest 
E T are labelled "1" and "2". 
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for 4 > 5 QeV and -2.5 < 77W < 0, i = 1, 2, (x 7 ) ~ 0.25, is just the region where pointlike gluons 
dominate D e ^(x, P 2 , M 2 ). This makes jet production in the region P 2 > 1 a promising place for 
identification of nontrivial aspects of PDF of virtual photons. 



4 PDF of the virtual photon in NLO QCD calculations 

Most of the existing information on interactions on virtual photons comes from the measurements 
of F 2 7 at PETRA |1| and LEP @ and jet production in ep collisions at HERA In [gl| 

data on dijet production in the region of virtualities 1 < P 2 < 80 GeV 2 , and for jet transverse 
energies E™ > 5 GeV have been analyzed within the framework of effective PDF defined in (p5|). 
This analysis shows that in the kinematical range 1 GeV 2 <P 2 <^E 2 the data agree reasonably 
with the expectations based on SaS parameterizations of PDF of the virtual photon. The same 
data may, however, be also analyzed using the NLO parton level Monte-Carlo programs that do 
not introduce the concept of PDF of the virtual photon. Nevertheless, so long as P 2 <C M 2 ~ Ej,, 
the pointlike parts of PDF incorporate numerically important effects of a part of higher order 
corrections, namely those coming from collinear emission of partons in Fig. |]. This makes the 
concept of PDF very useful phenomenologically even for the virtual photon. To illustrate this 
point we shall now discuss dijet cross-sections calculated by means of JETVIP [25|, currently the 



only NLO parton level MC program that includes both the direct and resolved photon contributions 
and which can thus be used to investigate the importance of the latter. 

4.1 Structure of JETVIP 

All the above mentioned parton level NLO MC programs contain the same full set of partonic cross- 
sections for the direct photon contribution up the order aa 2 . Examples of such diagrams P] are in 
Fig. |l2| a (a(aa s ) tree diagram) and Fig. |l2^ b (a(aa 2 ) tree diagram). Moreover all these programs 
contain also one-loop corrections to a(aa s ) tree diagrams. They differ mainly in the technique 
used to regularize mass singularities: MEPJET and JETVIP employ the slicing method whereas 
DISENT and DISASTER use the subtraction method. Numerical comparison of JETVIP and the 
other codes can be found in |f||7). To go one order of a s higher and perform complete calculation 



of the direct photon contributions up to order aa s would require evaluating tree diagrams like those 



in Fig. 12f,k, as well as one-loop corrections to diagrams like in Fig. 12 b and two-loop corrections 



to diagrams like in Fig. 12a. So far, such calculations are not available. 

In addition to complete NLO direct photon contribution JETVIP includes also the resolved 
photon one. Once the concept of virtual photon structure is introduced, part of the direct photon 
contribution, namely the splitting term (|2"3|), and in higher orders also further terms in (]T8|), is 
subtracted from the direct contribution (which for the virtual photon is nonsingular) and included 
in PDF appearing in the resolved photon contribution. To avoid confusion we shall henceforth use 
the term "direct unsubtracted" (DIR uns ) to denote NLO direct photon contributions before this 
subtraction and reserve the term "direct" for the results after it. In this terminology the complete 
calculations is then given by the sum of direct and resolved parts and denoted DIR+RES. 



10 Three such programs do exist, DISENT, Q, MEPJET, |23| and DISASTER @. 

11 In this subsection the various terms considered will be characterized by the powers of a and a s that appear in hard 
scattering cross-sections. In the corresponding Feynman diagrams of Fig. ^ these powers are given by the number 
of electromagnetic and strong vertices. Writing a(a? a 3 ) will thus mean parton level cross-sections proportional to 
cr'aj, not terms up to this order! For the latter we shall employ the standard symbol O(oyas). Because PDF of the 
photon are proportional to a, their convolutions in the resolved channel with partonic cross-sections cr(aj') are of the 
same order as partonic cross-sections a-(aa^) in the direct channel. For approximations taking into account the first 
two or three powers of a s , in either direct or resolved channel, the denomination NLO and NNLO will be used. 
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At the order a 2 the addition of resolved photon contribution means including diagrams like 
those in Fig. |l^c-e, which involve convolutions of PDF from both proton and photon sides with 
c(ct 2 ) tree partonic cross-sections. For a complete 0(a 2 ) calculation this is all that has to be added 
to the 0(aa 2 ) partonic cross-sections in direct photon channel. However, for reasons discussed in 
detail in the next subsection, JETVIP includes also NLO resolved contributions, which involve 
convolutions of PDF with complete cr{ot^.) partonic cross-sections (examples of relevant diagrams 
are in Figs. |l^g-j). This might seem inconsistent as no corresponding a(aaf) direct photon terms 
are included. Nevertheless, this procedure makes sense precisely because of a clear physical meaning 
of PDF of the virtual photon! Numerically, the inclusion of the cr(af) resolved terms turns out to 
be very important and in certain parts of the phase space leads to large increase of JETVIP results 
compared to those of DISENT, MEPJET or DISASTER. 



4.2 Factorization mechanism in 7p interactions 

The main argument for adding <r(a^) partonic cross-sections in the resolved channel to 0{aa 2 s ) 
ones in the direct and 0(a 2 ) ones in the resolved channels, is based on specific way factorization 
mechanism works for processes involving photons in the initial state. First, however, let us recall 
how the factorization works in hadronic collisions. Jet cross-sections start as convolutions of &(a 2 ) 
partonic cross-sections with PDF of beam and target hadrons. The factorization scale dependence 
of these PDF is cancelled by explicit factorization scale dependence of higher order partonic cross- 
sections. This cancellation is exact provided all orders of perturbation theory [^] are taken into 
account, but only partial in any finite order approximation, like the NLO one used in analyses of 
jet production at FERMILAB. In hadronic collisions the inclusion of cr(a^) partonic cross-sections is 
thus vital for compensation of the factorization scale dependence of PDF in convolutions with <r(a 2 ) 
partonic cross-sections. The residual factorization scale dependence of such NLO approximations 
is formally of the order and thus one order of a s higher than the terms included in the NLO 
approximation . 

In 7p collisions (whether of real or virtual photon) the situation is different due to the presence 
of inhomogeneous terms k q ,kc in the evolution equations (|3]-||). For quark distribution functions 
(nonsinglet as well as singlet) the leading term {a/2ir)k q ^ in the expansion of k q is independent 
of a s and consequently part of the factorization scale dependence of the cr(a 2 ) resolved photon 
contribution is compensated at the same order a 2 . In further discussions we shall distinguish two 
factorization scales: one (M 7 ) for the photon, and the other {M p ) for the proton. The content 
of the evolution equations (H|-||) is represented graphically in Fig. |l2|, which shows examples of 
diagrams, up to order , relevant for our discussion. Some of these diagrams are connected by 
solid or dashed arrows, representing graphically the effects of variation of factorization scales M 7 
and M p respectively. The vertical dashed arrows connect diagrams (with lower blobs representing 
PDF of the proton, denoted D p (x, M p )) that differ in partonic cross-sections by one order of a s , 
reflecting the fact that the terms on the r.h.s. of evolution equation for D p (x,M p ) start at the 
order a s . For instance, the a(aa s ) direct photon diagram in Fig. |l^ a is related by what we call 
Mp-factorization with a(aa 2 ) direct photon diagram in Fig. |l^b. Similarly, the <j{a 2 ) resolved 
photon diagram in Fig. |i~2|d is related by M p -factorization to <t(q^) resolved photon diagram in Fig. 



12h. In fact each diagram of the order is related by M„-factorization to two types of diagrams 



at order aj +1 , one with quark and the other with gluon coming from the proton blob. Note that 
Mp-factorization operates within either direct or resolved contributions separately, never relating 
one type of terms with the other. 



12 In perturbation expansions of partonic cross-sections as well as splitting functions. 



18 




G p (x,M p ) 



Figure 12: Examples of diagrams related by factorization mechanism. Only powers of a s in parton 
level cross-sections are counted. 
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For M 7 -factorization this cancellation mechanism has a new feature. Similarly as for hadrons, 



the cr(a s ) resolved photon diagram in Fig. 12c is related by what we call homogeneous M 7 



factorization to the cr(af) resolved photon diagrams in Fig. [l2]g,h, and similar relation holds also 
between the diagram in Figs. |l2| e and 12 j. However, the inhomogeneous term in evolution equation 
for PDF of the photon implies additional relation (which we call inhomogeneous M 7 -factorization) 
between direct and resolved photon diagrams at the same order of a s , represented by horizontal 
solid arrows. For instance, the LO resolved contribution coming from diagram in Fig. 12] c is 
related not only by homogeneous M 7 -factorization to the cr(ag) resolved photon diagrams L|g-h, 
but also by inhomogeneous M 7 -factorization to the a(aa 2 ) direct photon diagram in Fig. |l2|b. 
Similarly, the cr(a^) resolved photon diagram p^ g is related by homogeneous M 7 -factorization to 
o~{a^) resolved terms (not shown) and by inhomogeneous M 7 -factorization to the a(aa^) direct 



diagram in Fig. |l2| f. 

The argument for adding the cr(a^) resolved photon diagrams to the o~{a 2 ) resolved photon 
ones relies on the fact that within the resulting set of NLO resolved contributions the homogeneous 
M 7 -factorization operates in the same way as the M p -factorization does within the set of NLO 
direct ones! However, contrary to the hadronic case, the o~(a? s ) parton level cross-sections do 
not constitute after convolutions with photonic PDF a complete set of aa^ contributions, the rest 
coming from the cr(aa^) direct ones. As these a(aa^) direct photon contributions have not yet been 
calculated, the set of contributions included in JETVIP does not constitute the complete 0{ao? s ) 
calculation. The lacking terms, like that coming from the diagram 12f, would provide cancellation 
mechanism at the order aa^ with respect to the inhomogeneous M 7 -factorization. In the absence 
of a(aa^) direct photon calculations, we thus have two options: 



• To stay within the framework of complete 0(a 2 ) calculations, including the LO resolved and 
NLO direct contributions, but with no mechanism for the cancellation of the dependence of 
PDF of the virtual photon on the factorization scale M 7 . 

• To add to the previous framework the cr(a^) resolved photon contribution, which provide the 
necessary cancellation mechanism with respect to homogeneous M 7 -factorization, but do not 
represent a complete set of o"(af) contributions. 



In our view the second strategy, adopted in JETVIP, is more appropriate. In fact one can look at 
cr(a^) resolved photon terms as results of approximate evaluation of the so far uncalculated a(aa^) 
direct photon diagrams in the collinear kinematics. For instance, so long as P 2 <C M 2 taking into 
account only the pointlike part of D 7 (x, P 2 , M 2 ) in the upper blob of Fig. |l2] g should be a good 
approximation of the contribution of the direct photon diagram in Fig. f[2] f. There are of course 
o~(aa^) direct photon diagrams that cannot be approximated in this way, but we are convinced 
that it makes sense to build phenomenology on this framework. 

For the o~{a 2 ) resolved terms the so far unknown a(aaf) direct photon contributions provide 
the first chance to generate pointlike gluons inside the photon: the gluon in the upper blob in 
resolved photon diagram in Fig. Il2fe must be radiated by a quark that comes from primary 7 — > qq 



splitting, for instance as shown in Fig. 12j. Note that to get the gluon in <r(a^) resolved photon 
contributions, for instance in Fig. 12h, within collinear limit of direct photon contributions would 
require evaluating diagrams, like that in Fig. 12m, at the order a(aaj,)\ In summary, although 
the pointlike parts of quark and gluon distribution functions of the virtual photon are in a sense 
included in higher order perturbative corrections and can therefore be considered as expressions of 
"interactions" rather than "structure" of the virtual photon, their uniqueness and phenomenological 
usefulness definitely warrant their introduction as well as name. 
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4.3 Theoretical uncertainties 



NLO calculations of jet cross-sections are affected by a number of ambiguities caused by the 
truncation of perturbation expansions as well as by uncertainties related to the input quantities 
and conversion of parton level quantities to hadron level observables. 



4.3.1 Choice of PDF 

We have taken CTEQ4M and SAS1D sets of PDF of the proton and photon respectively as our 
principal choice. Both of these sets treat quarks, including the heavy ones, as massless above their 
respective mass thresholds, as required by JETVIP, which uses LO and NLO matrix elements of 
massless partons. PDF of the proton are fairly well determined from global analyses of CTEQ 
and MRS groups and we have therefore estimated the residual uncertainty related to the choice of 
PDF of the proton by comparing the CTEQ4M results to those obtained with MRS(2R) set. The 
differences are very small, between 1% at r\ = —2.5 and 3.5% at r\ = 0, independently of P 2 . 

For the photon we took the Schuler-Sjostrand parameterizations for two reasons. First, they 
provide separate parameterizations of VDM and pointlike components of all PDF, which is crucial 
for physically transparent interpretation of JETVIP results. Secondly, they represent the only set 
of photonic PDF with physically well motivated virtuality dependence, which is compatible with 
the way JETVIP treats heavy quarks. The GRS sets |14,15|, the only other parameterizations with 



built-in virtuality dependence, are incompatible with JETVIP because they treat c and b quarks 
as massive and, consequently, require calculating their contribution to physical quantities via the 
boson gluon or gluon gluon fusion involving exact massive matrix elements. 



4.3.2 The choice of n/ 

JETVIP, as well as other NLO parton calculations of jet cross-sections, works with a fixed number 
rif of massless quarks, that must be chosen accordingly. This is not a simple task, as the number of 
quarks that can be considered effectively massless depends on kinematical variables characterizing 
the hardness of the collision. Consequently, the optimal choice of nj may not be unique for the 
whole kinematical region under consideration. The usual procedure is to run such programs for 
two (or more) relevant values of and use the ensuing difference as an estimate of theoretical 
uncertainty related to the approximate treatment of heavy quark contributions. 

The number nj enters NLO calculations in three places: implicitly in and PDF and 

explicitly in LO and NLO parton level cross-sections. In our selected region of phase space the 
appropriate value of nj lies somewhere between rif = 4 and rif = 5, with the latter value represent- 
ing the upper bound on the results (so far unavailable) that would take the fe-quark mass effects 
properly into account. We have therefore run JETVIP for both n/ = 4 and rif = 5 and compared 
their results. They differ, not suprisingly, very little 0. Explicit calculations give at most 10% 
in the direct channel and 5% in the resolve done. All the results presented below correspond to 
n f = 5. 

4.3.3 Factorization and renormalization scale dependence 

As mentioned in the previous subsection, proton and photon are associated in principle with 
different factorization scales M p and M 7 , but we followed the standard practice of assuming 

As the contribution of the 6-quark in both direct and resolved channels is proportional to its charge, we except 
it to amount to about ej/(e^ + + + e^) = 0.1 of the sum of light quark ones. 
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M = M p = M-y and set M = kEj, . The factorization scale dependence was quantified by 
performing the calculations for k = 0.5,1 and n = 2. 

The dependence of finite order perturbative calculations on the renormalization scale \x is in 
principle a separate ambiguity, but we have again followed the common practice of identifying 
these two scales /J, = M. To reflect this identification, we shall in the following use the term "scale 
dependence" to describe the dependence on this common scale. 



4.3.4 Hadronization corrections 

JETVIP, as a well as other NLO codes evaluate jet cross-sections at the parton level. For a mean- 
ingful comparison with experimental data they must therefore be corrected for effects describing 
the conversion of partons to observable hadrons. These so called hadronization corrections are not 
simple to define, but adopting the definition used by experimentalists [28] we have found that they 



depended sensitively and in a correlated manner on transverse energies and pseudorapidities of 
jets. In order to avoid regions of phase space where they become large we imposed on both jets the 
condition —2.5 < 77 W < 0. In this region hadronization corrections are flat in r\ and do not exceed 
10%, whereas for rj = —2.5 they steeply rise with decreasing ij. This by itself would not require 
excluding this region, the problem is that in this region hadronization corrections become also very 
much model dependent and therefore impossible to estimate reliably. Detailed analysis of various 
aspects of estimating these hadronization , with particular emphasize on their implication for jet 



production at HERA, is contained in 2£]. 



4.3.5 Limitations of JETVIP calculations 

Despite its undisputable advantage over the calculations that do not introduce the concept of virtual 
photon structure, also JETVIP has a drawback because it does not represent a complete NLO QCD 
calculation of jets cross-sections. This is true in the conventional approach to photonic interactions, 
and even more in the reformulation suggested by one of us in Q . In the conventional approach the 
incompletness is related to the fact that there is no NLO parameterization of PDF of the virtual 
photon compatible with JETVIP treatment of heavy quarks. Note that in the standard approach 
the inclusion of cr(a^) partonic cross-sections in the resolved photon channel is justified by the claim 
that their convolution with photonic PDF are of the same order aa 2 as the direct photon ones. 
In the reformulation |l| this incompletness has deeper causes. It reflects the lack of appropriate 
input PDF but also the fact that a complete NLO approximation requires the inclusion of direct 
photon contribution of the order aa^, which is so far not available. Nevertheless, we reiterate that 
it makes sense to build phenomenology upon the current JETVIP framework and the concept of 
PDF of the virtual photon is just the necessary tool for accomplishing it. 



4.4 Dijet production at HERA 

We shall now discuss the main features of dijet cross-sections calculated by means of JETVIP. To 
make our conclusions potentially relevant for ongoing analyses of HERA data we have chosen the 
following kinematical region 

Eft > E C T + A, Eft > E%, E^ = 5 GeV, A = 2 GeV 

-2.5 < r/ W < 0, i = 1,2 

in four windows of photon virtuality 

1.4 < P 2 < 2.4 GeV 2 ; 2.4 < P 2 < 4.4 GeV 2 ; 4.4 < P 2 < 10 GeV 2 ; 10 < P 2 < 25 GeV 2 
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The cuts on Et were chosen in such a way that throughout the region P 2 <C E T , thereby ensuring 
that the virtual photon lives long enough for its "structure" to develop before the hard scattering 
takes place. The asymmetric cut option is appropriate for our decision to plot the sums of Et and 
rj distributions of the jets with highest and second highest Et- The choice of A = 2 GeV, based on 
a detailed investigation |H| of the dependence of the integral over the selected region on A, avoids 
the region where this dependence possesses unphysical features. 

In our analysis jets are defined by means of the cone algorithm. At NLO parton level all 
jet algorithms are essentially equivalent to the cone one, supplemented with the parameter R SC p, 



introduced in [30 in order to bridge the gap between the application of the cone algorithm to 
NLO parton level calculations and to hadronic systems (from data or MC), where one encounters 
ambiguities related to seed selection and jet merging. In a general cone algorithm two objects 
(partons, hadrons or calorimetric cells) belong to a jet if they are within the distance R from the 
jet center. Their relative distance satisfies, however, a weaker condition 



AR 



„ - \J (Ar/ij) 2 + (A^) 2 < J^ E f \ R - ( 46 ) 
v max(£/Ti , hrr. ) 



The parameter i? S ep governs the maximal distance between two partons within a single jet, i.e. two 
partons form a jet only if their relative distance ARij satisfies the condition 



ARij < min 



--TL, -Tt S ep 



max(E Tl , E Ti ) ' 



(47) 



The question which value of -R se p to choose for the comparison of NLO parton level calculations 
with the results of the cone algorithm at the hadron level is nontrivial and we shall therefore present 
NLO results for both extreme choices R se p = R and -R se p = 2i?. To define momenta of jets JETVIP 
uses the standard .Ey-weighting recombination procedure, which leads to massless jets. 

4.5 Results 

To asses phenomenological importance of the concept of PDF of virtual photons we now compare 
JETVIP results obtained in the DIR uns mode, where this concept is not introduced at all, with 
those of the DIR+RES one, in which the contribution of the resolved photon is added to the 
subtracted direct one. The difference between these two results measures the nontrivial aspects of 
PDF of the virtual photon. 



We start by plotting in Fig. 13 the distributions da/drj and da/dET in the first window 
1.4 < P 2 < 2.4 GeV 2 . All curves correspond to R SC p = 2. The difference between the solid and 
dashed curves is significant in the whole range of 77, but becomes truly large close to the upper edge 
rj = 0, where the DIR+RES results exceed the DIR uns ones by a factor of about 3! In da/dET 
distributions this difference comes predominantly from the region of Et close to the lower cut-off 



E T + A = 7 GeV. Fig. |13| also shows that the scale dependence is nonnegligible for both DIR uns 
and DIR+RES results, but does not invalidate the main conclusion drawn from this comparison. 
Interestingly, the scale dependence is weaker for the DIR+RES results than for the DIR uns ones. 
Fig. [14] documents that the above results depend only very weakly on R sep . 



To track down the origins of the observed large differences between DIR+RES and DIR 



uns 



results we compare in Fig. [1^ the DIR+RES and DIR uns results to the subtracted direct ones 
(denoted DIR). The difference between the DIR+RES and DIR curves, giving the resolved photon 
contribution da res /drj, is further split into the following contributions: 

• VDM part of photonic PDF convoluted with complete NLO parton level cross-sections (de- 
noted NLO VDM). 
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Figure 13: Scale dependence of the distributions da/drj and da/dEx at the NLO. All curves 
correspond to R sep = 2R. 




Figure 14: i? sep dependence of da/drj and da/dEr distributions. All curves correspond to fi = Et- 
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Figure 15: Comparison of DIR+RES, DIR uns and DIR results for da/dr] (left plots) and individual 
contributions to da Tes /dr], described in the text (right plots). The thin solid curve corresponds to 
convolution of the splitting term (23) with a v g es (o?) parton level cross-sections. 
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Figure 16: Fractional contributions to d<7 res /d?7. Upper and lower dotted (dashed-dotted) curves 
correspond to pointlike quarks (gluons) convoluted with 0{o? s ) and 0{a 2 s ) partonic cross-sections. 
The dashed curve in a) corresponds to the NLO VDM contribution. The solid curves denote the 
ratio DIR uns /DIR+RES . 
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Figure 17: Nontriviality fractions R3 and R4 as functions of rj and P 2 



• Pointlike quarks and gluons convoluted with a ves (a 2 ) and a ies (a 3 ) parton level cross-sections 

displayed separately in the upper right plot of Fig. The full NLO resolved photon contribution 
is given as the sum 



NLO VDM + ai es (a 2 ) <g> q PL + a r ^ s (a 2 ) <g> G PL + a* es (a 3 ) <g> q PL + cr£f(a s 3 ) G 



PL 



(48) 



Fractional contributions of LO and NLO terms to <r res coming separately from pointlike quarks 
and gluons are plotted in Fig. |i"6| a as functions of r\. Several conclusions can be drawn from Figs. 
T5HT1 



• The contribution of the VDM part of photonic PDF is very small and perceptible only close 
to r\ = 0. Integrally it amounts to about 3%. Using SaS2D parameterizations would roughly 
double this number. 

• The inclusion of al es (a 3 ) parton level cross-sections in the resolved photon channel is numer- 
ically very important throughout the range —2.5 < rj < 0. Interestingly, the o"[ es (a^) results 
are close, particularly for the pointlike quarks, to the ^^(a. 2 .) ones. 

• At both a 2 and a 3 orders pointlike quarks dominate da res /dr] at large negative rj, whereas as 
rj — > the fraction of d<T rcs /dry coming from pointlike gluons increases towards 40% at r/ = 0. 

We emphasize that pointlike gluons carry nontrivial information already in convolutions with 
a ies (a 2 ) partonic cross-sections because in unsubtracted direct calculations such contributions 
would appear first at the order aa 3 , [^]. The convolution of the dominant part of pointlike quarks p| 
with a ies (a 3 ) partonic cross-sections would be included in direct unsubtracted calculations starting 
also at the order aa 3 , whereas for pointlike gluons this would require evaluating the unsubtracted 

14 For instance, the resolved photon diagram in Fig. p^ e would come as part of the of evaluating the unsubtracted 
direct diagram in Fig. [l2|k, 

lj That is, the one given by the splitting term (|23|). 



26 



direct terms of even higher order aaj.1 For instance, the contribution of diagram in Fig. 12 g would 



be included in the contribution of diagram in Fig. 12 f. Similarly, the results of diagram in Fig. 



12h would come as part of the results of evaluating the diagram in Fig. [12;m. 

In JETVIP the nontrivial aspects of taking into account the aj es (ag) resolved photon contribu- 
tions can be characterized [^] by the "nontriviality fractions" R3 and R4 

^ g PL ®<>f) + G PL ®ag s (a 2 ) _ G PL ^^(q3) 

3 <j rcs ' <7 res 

which quantify the fractions of <r res that are not included in NLO unsubtracted direct calculations. 
These fractions are plotted as functions of r\ and P 2 in Fig. [l7|. Note that at 7] = almost 70% of 
cr res comes from these origins. This fraction rises even further in the region rj > 0, which, however, 
is experimentally difficult to access. 

So far we have discussed the situation in the first window of photon virtuality, i.e. for 1.4 < 
P 2 < 2.4 GeV 2 . As P 2 increases the patterns of scale and i? S ep dependencies change very little. 



On the other hand, the fractions plotted in Fig. 16 and 17 vary noticeably: 



• The DIR uns contributions represent an increasing fractions of the DIR+RES results. 

• The relative contribution of pointlike gluons with respect to pointlike quarks decreases. 

• The nontriviality factor R4 (which comes entirely from pointlike gluons) decreases, whereas 
-R3, which is dominated by pointlike quarks and flat in 77, is almost independent of P 2 . 

All these features of JETVIP results reflect the fundamental fact that as P 2 rises towards the 
factorizations scale M 2 ~ E 2 -, the higher order effects incorporated in pointlike parts of photonic 
PDF vanish and consequently the unsubtracted direct results approach the DIR+RES ones. The 
crucial point is that for pointlike quarks and gluons this approach is governed by the ratio P 2 /M 2 
appearing in the multiplicative factor (1 — P 2 /M 2 ). The nontrivial effects included in PDF of 
the virtual photon will thus persist for arbitrarily large P 2 , provided we stay in the region where 
P 2 < M 2 . Moreover, they are so large, that they should be visible already in existing HERA 
data. Provided the basic ideas behind the Schuler-Sjostrand parameterizations of PDF of the 
virtual photon are correct, our analysis shows that the calculations that do not introduce the 
concept of virtual photon structure should significantly undershoot the available HERA data on 
dijet production in the kinematical region 1 < P 2 < E 2 , E T >5 GeV, -2.5 < 77 < 0, in particular 
for 77 ~ 0. Published as well preliminary data discussed in [29,31] support this conjecture. 



5 Summary and conclusions 

We have analyzed the physical content of parton distribution functions of the virtual photon within 
the framework formulated by Schuler and Sjdstrand, which provides physically motivated separation 
of quark and gluon distribution functions into their hadronic (VDM) and pointlike parts. We have 
shown that the inherent ambiguity of this separation, numerically large for the real photon, becomes 
phenomenologically largely irrelevant for virtual photons with P 2 > 2 — 3 GeV 2 . In this region 
quark and gluon distribution functions of the virtual photon are dominated by their (reasonably 
unique) pointlike parts, which have clear physical origins. We have analyzed the nontrivial aspects 
of these pointlike distribution functions and, in particular, pointed out the role of pointlike gluons 
in leading order calculations of jet cross-section at HERA. 

16 Disregarding the VDM part of resolved contribution which is tiny in our region of photon virtualities. 
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The conclusions made within the framework of LO QCD have been confirmed, and in a sense 
even strengthened, in our analysis of NLO parton level calculations using JETVIP. We have found 
a significant difference between JETVIP results in approaches with and without the concept of 
virtual photon structure. While for the real photon analogous difference is in part ascribed to 
the VDM part of photonic PDF, for moderately virtual photons it comes almost entirely from 
the pointlike parts of quark and gluon distribution functions. Although their contributions are in 
principle contained in higher order calculations which do not use the concept of PDF, in practice 
this would require calculating at least a(ao? s ) and a(aa^) unsubtr acted direct contributions. In 
the absence of such calculations the concept of PDF of the virtual photon is therefore very useful 
phenomenologically and, indeed, indispensable for satisfactory description of existing data. 

Acknowledgment: We are grateful to J. Cvach, J. Field, Ch. Friberg, G. Kramer, B. Potter, 
I. Schienbein and A. Valkarova for interesting discussions concerning the structure and interactions 
of virtual photons and to B. Potter for help in running JETVIP. This work was supported in part 
by the Grant Agency of the Academy od Sciences of the Czech Republic under grant No. A1010821. 

References 



[i] 

[2] 
[3] 
[4] 
[5] 

[6] 

[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 



J. Chyla: jhep-ph/99114T3| 

W. Bardeen, A. Buras, D. Duke and T. Muta, Phys. Rev D18, 3998 (1978) 
J. Chyla, Z. Phys. C43, 431 (1989) 

G. Curci, W. Furmanski, R. Petronzio, Nucl. Phys. B175, 27 (1980) 

P. Aurenche, R. Baier, A. Douiri, M. Fontannaz, D. Schiff, Nucl. Phys. B286, 553 (1987) 
P. Aurenche, R. Baier, M. Fontannaz, D. Schiff, Nucl. Phys. B296, 661 (1987) 

J. Chyla, M. Tasevsky, in Proceedings PHOTON '99, Freiburg in Breisgau, May 1999, ed. S. 



Soeldner-Rembold, Nucl. Phys. Proc. Supp., in press, [hep-ph / 9906552 



G. Schuler, T. Sjostrand: Z. Phys. C68 (1995), 607 
G. Schuler, T. Sjostrand: Phys. Lett. B376 (1996), 193 

J.H. Field, F. Kapusta, L. Poggioli, Phys. Lett. B181, 362 (1986) 
J.H. Field, F. Kapusta, L. Poggioli, Z. Phys. C36, 121 (1987) 
F. Kapusta, Z. Phys. C42, 225 (1989) 

E. Witten: Nucl. Phys. 120 (1977). 189 

W. Bardeen, A. Buras, Phys. Rev. D20, 166 (1979) 

J. Chyla: Phys. Lett. B320 (1994), 186 

J. Chyla: [hep-ph/ 98 11455 



A. Gorski, B.L. Ioffe, A. Yu. Khodjamirian, A. Oganesian, Z. Phys. C 44, 523 (1989) 

M. Gliick, E. Reya, M. Stratmann, Phys. Rev. D51 (1995), 3220 
M. Gliick, E. Reya, M. Stratmann, Phys. Rev. D54 (1996), 5515 



[15] M. Gliick, E. Reya, I. Schienbein, Phys. Rev. D60 (1999), 054019 



28 



[16] 
[17] 

[18] 
[19] 

[20] 
[21] 
[22] 
[23] 
[24] 
[25] 

[26] 

[27] 
[28] 

[29] 
[30] 
[31] 



J. Chyla, M. Tasevsky, in preparation 

J. Chyla, M. Tasevsky, in Proceedings of Workshop MC generators for HERA Physics, Ham- 
burg 1999, p. 239, |hep-ph/ 9905444 



Ch. Berger et al. (PLUTO Collab.), Phys. Lett. B142 (1984), 119 

F. Erne, in Proceedings PHOTON '99, Freiburg in Breisgau, May 1999, ed. S. Soeldner- 
Rembold, Nucl. Phys. Proc. Supp., in press, L3 Note 2403, submitted to Europhysics Confer- 
ence High Energy Physics, Tampere, August 1999 

C. Adloffet al. (HI Collab.), Phys. Lett. B415 (1997), 418 

C. Adloff et al. (HI Collab.), Eur. Phys. J. C, in press 

S. Catani, M. Seymour, Nucl. Phys. B458 (1997), 291, Erratum-ibid. B510 (1997), 503 
E. Mirkes, TTP-97-39, |hep-ph/97li22l 



D. Graudenz, |hep-ph/9710244 



B. Potter, Comp. Phys. Comm. 119 (1999), 45, |hep-ph/9806~437 
G. Kramer, B. Potter, Eur. Phys. J. C5 (1998), 665, |hep-ph/9804352 



C. Duprel, Th. Hadig, N. Knauer, M. Wobisch, in Proceedings of Workshop MC generators for 



HERA Physics, Hamburg 1999, p. 142, |hep-ph/991044g 



B. Potter, |hep-ph/9911221 



M. Wobisch, in Proceedings of Workshop MC generators for HERA Physics, Hamburg 1999, 



p. 239, |hep-ph/9905444 

M. Tasevsky, PhD Thesis, unpublished 

S.D.Ellis, Z. Kunszt, D.E. Soper, Phys. Rev. Lett. 69 (1992), 3615 

J. Cvach, in Proceedings DIS99, Nucl. Phys. B (Proc. Suppl.) 79 (1999), 501 



29 



